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O ' 

' We propose a practical method for analyzing stability of Q-balls for the whole parameter 

space, which includes the intermediate region between the thin-wall limit and thick-wall 
limit as well as Q-bubbles (Q-balls in false vacuum), using catastrophe theory. We apply 
our method to the two concrete models, V3 = rn 2 (f> 2 /2 — (iff) 3 + Xcj> 4 and V& = m 2 (j> 2 /2 — 
A0 4 + <t> 6 /M 2 . We find that V3 and Va Models fall into fold catastrophe and cusp catastrophe, 
\ respectively, and their stability structures are quite different from each other. 
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Q-balls, 1 ^ a kind of non-topological solitons, 2 ) appear in a large family of field 
theories with global U(l) (or more) symmetry, and could play an important role in 
cosmology. For example, the Minimal Supersymmetric Standard Model may con- 
tain baryonic Q-balls, which could be responsible for baryon asymmetry 3 ) and dark 
matter. 4 ) 

psj . The stability of Q-balls has been studied in the literature. Coleman argued that 

Q-balls are absolutely stable if the charge Q is sufficiently large, using the thin-wall 
approximation. 1 ) Kusenko showed that Q-balls with small Q are also stable for the 
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Csj! W) = ^-0 2 - W> 3 + A^ 4 with m 2 , /i, A>0, (1) 

t-h ; 

using the thick-wall approximation. 5 ) Here the thick-wall limit is defined by the 
limit of to 2 — ► m 2 , where to is the angular velocity of phase rotation. Multamaki 
and Vilja found that in the thick-wall limit the stability depends on the form of the 
^ ■ potential. 6 ) Paccetti Correia and Schmidt showed a useful theorem which applies to 

any equilibrium Q-balls: 7 ) their stability is determined by the sign of {uj / Q)dQ / doo . 

It is usually assumed that the potential has an absolute minimum at (j) = 0. If 
V(0) is a local minimum and the absolute minimum is located at <f> 7^ 0, true vacuum 
bubbles may appear. 8 ) If Q = 0, vacuum bubbles are unstable: either expanding or 
contracting. Kusenko 9 ) and Paccetti Correia and Schmidt 7 ) showed, however, that 
there are stable bubbles if Q 7^ 0. They called those solutions "Q-balls in the false 
vacuum" . Hereafter we simply call them "Q-bubbles" . 

The standard method for analyzing stability is to take the second variation 
of the total energy (given by Eq.Q below) and evaluate its sign. However, this 
calculation can be executed analytically only for some limited cases; in general the 
eigenvalue problem should be solved numerically, as Axenides et al. did. 10 ) In this 
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paper, we propose an easy and practical method for analyzing stability with the help 
of catastrophe theory. The basic idea of catastrophe theory is described in Appendix. 
As we shall show below, once we find behavior variable(s), control parameters) and 
a potential in the Q-ball system, it is easy to understand the stability structure of 
Q-balls for the whole parameter space including the intermediate region between the 
thin-wall limit and thick-wall limit as well as Q-bubbles. 

Consider an S0(2)-symmetric scalar field, whose action is given by 



d 4 x 



with <f> = J $ + 4%. (2 



We consider spherically symmetric configurations of the field. Assuming homoge- 
neous phase rotation, 

(</>i> 4>2) = 4>(r)(cos ut, smut), (3) 

the field equation becomes 

^ = _2#_ w 2 + ^ (4) 
dr 2 r dr d4> 

This is equivalent to the field equation for a single static scalar field with the potential 
V W = V — u 2 </> 2 /2. Due to the symmetry there is a conserved charge, 

Q = J d A x{(t} 1 d t (t)2 ~ 02<9t0i) = ojI, where I = Jd 3 x<p 2 . (5) 

Monotonically decreasing solutions (j)(r) with the boundary conditions, 

^(0)=0, 0(oo) = 0, (6) 
dr 

exist if min(V^) < V(0) and d 2 V u) /d(f) 2 (0) > 0, which is equivalent to 



UJ mm< uj2 <m 2 with u4 in = min — v , m 2 = — v (0) , (7) 



mm 



'2V\ 2 = (fV 
' mm ^ ^ ^ "" " mm ~ V ~W )' m = ~d~4 2 

where we have put V(0) = without loss of generality. The two limits uj 2 — > uj 2 
and uj 2 — > m 2 correspond to the thin-wall limit and the thick-wall limit, respectively. 
The condition a ^ in < m 2 is not so restrictive because it is satisfied if the potential 
has the form, 

2 

V = ^-cp 2 - \<p n + 0{(f> n+1 ) with m 2 > 0, A > 0, n > 3 . (8) 
The total energy of the system for equilibrium solutions is given by 

Q 2 , /\ 3 „fl (d^ 1 



Note that the variation of E under fixed Q, 5E/54>\q = 0, reproduces the field 
equation (jH). 
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Let us discuss how we apply catastrophe theory to the present Q-ball system. 
Catastrophe theory is briefly described in Appendix. An essential point is to choose 
behavior variable(s), control parameters) and a potential in the Q-ball system ap- 
propriately. For a given potential V(<j>) and charge Q, we consider a one-parameter 
family of perturbed field configurations (j>u(r) near the equilibrium solution 0(r). 
The one-parameter family is chosen to satisfy I[(j>w] = Q/u>. Then the energy is 
regarded as a function of u, E(uS) = E\<j)J\. 

Because dE/duj = {5E / 5(f)^)d^) w / a\o = when W is an equilibrium solution, u 
may be regarded as a behavior variable and E as the potential. On the other hand, 
the charge Q and the model parameter(s) of V(4>) can be given by hand, and therefore 
should be regarded as control parameters. We denote the model parameter(s) by Pi 
(i = 1, 2, • • • ). Then we analyze the stability of Q-balls as follows. 

• Solve the field equation with the boundary condition ([6]) numerically to 
obtain equilibrium solutions 0(r) for various values of to and model parameter(s) 
Pi- 

• Calculate Q by (|5|) for each solution to obtain the equilibrium space M = 
{(to, Pi, Q)}. We denote the equation that determines M by f(tu, Pi, Q) = 0. 

• Find folding points where dPi/duj = or dQ/du = in M, which are identical 
to the stability-change points, U = {(to, Pi,Q) \ df/dto = 0, / = 0}. 

• Calculate the energy E by ([9]) for equilibrium solutions around a certain point 
in 27 to find whether the point is a local maximum or a local minimum. Then 
we find the stability structure for the whole M. 

Now, using the method devised above, we investigate the stability of equilibrium 
Q-balls. Because it was shown 7 * 1 that in the thick-wall limit Q-balls are stable if 
n < 10/3 for the potential ([8]) and unstable otherwise, we consider two typical 
models. One is given by (pQ), which we call V3 Model, and the other is given by 

m 2 df 3 
Vt{4>) = —4? - A0 4 + ^ with m 2 , A, M 2 > , (10) 

which we call V4 Model. For V3 Model, rescaling the quantities as 

T- + ~_ 1- X A t>_ a3 t/ ~ _ 

t = —t, r = —r, (f> = -0, V 3 = —TV3, m = m, to = lu, (11) 

VA VA M M A* A- 4 

the field equation (jU), the potential ([T]), the charge © and the energy ([9]) are rewrit- 
ten as 

~ ~ -* 3 

d 2 d> 2dd> 9 - dVz ~ rh 2 79 T o ~ A ~ A 2 ~ ^ N 

W = ~d-^* + -i' V 3 = ^ 2 -* 3 + * 4 . Q = AQ. (12) 

Similarly, for V4 Model, rescaling the quantities as 

t = XMt, f = XMr, 4)=—{L- , V A =-^-r, m = -^-, w = (13) 

v^M A 3 M 4 ' AM AM v ; 

the field equation (jlj), the potential (fTUj) . the charge ([5]) and the energy ([9]) are 
rewritten as 



— = --— -w + -y, V 4 = — 0-0+0°, E = —, Q = —. 
dr z r dr dd 2 MX 
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In both models the system is regarded as a mechanical system with the behavior 
variable Co, the control parameters m 2 and Q, and the potential E(Q; m 2 , Q). Because 
^min = r "' 2 — 1/2; the existing condition ([7]) reduces to 

0<m 2 -Co 2 <^. (15) 

The thin- wall and thick- wall limits correspond to m 2 — Co 2 — > 1/2 and m 2 — u) 2 — > 0, 
respectively. The condition for ordinary Q-balls, > 0, reduces to to 2 > 1/2, 

while that for Q-bubbles, u; 2 ^ < 0, to m 2 < 1/2. 

Figures [T]and [2] show the structures of the equilibrium spaces, M = {(Co, m 2 , Q)}, 
and their catastrophe map, x(M), into the control planes, C = {(m 2 ,Q)}> for V3 
and V4 Models, respectively. We only show the results for Co > 0; the sign trans- 
formation Co — > — tD changes nothing but Q — ► — Q- The dash-dotted lines in M 
denote stability-change points £. Because the equilibrium space alone does not tell 
us which lines, solid or dashed, represent stable solutions, we evaluate the energy E 
for several equilibrium solutions, as shown in Figs. [3] and HI When there are double 
or triple values of E for a given set of the control parameters (m 2 , Q), by energetics 
the solution with the lowest value of E should be stable and the others should be 
unstable. In Figs. [3]and[U we also give a sketch of the potential E(u; to 2 , Q) near the 
equilibrium solutions. Once the stability for a given set of the parameters (m 2 , Q) is 
found, the stability for all the sets of parameters which may be reached continuously 
from that set without crossing £ is the same. We therefore conclude that, in Figs. 
[Hand [2] as well as in Figs. [3]and[U solid and dashed lines correspond to stable and 
unstable solutions, respectively. 

According to the configurations of x(^) m the control planes in Figs. [T] and 
[2j we find that V3 Model falls into fold catastrophe while V4 Model falls into cusp 
catastrophe. In the control planes, the numbers of stable and unstable solutions for 
each (rh 2 ,Q) are represented by N, S, U, SU and SUU (see the figure captions for 
their definitions). Thus we find the stability structures of the two models are very 
different from each other. They are found as follows. 
V 3 Model 

• to 2 > 1/2: All equilibrium solutions are stable. 

• to 2 < 1/2 (Q-bubbles): For each to 2 there is a maximum charge, Q max , 
above which equilibrium solutions do not exist. For Q < Q ma , x , stable and 
unstable solutions coexists. It is interesting to note that stable Q-bubbles 
exist no matter how small Q is. 

F 4 Model 

• to 2 > 1/2: For each to 2 there is a minimum charge, Q m i n , below which 
equilibrium solutions do not exist. For Q > Q m \ n , stable and unstable 
solutions coexists. 

• to 2 < 1/2 (Q-bubbles): For each to 2 there is a maximum charge, Q m & x , 
as well as a minimum charge, Qmin, where stable solutions do not exist if 
Q < Qmin or Q > Qmax- For Qmin < Q < Qmax, there are one stable and 
two unstable solutions. 
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Fig. 1. Structures of the equilibrium spaces, M = {(a), rh 2 , Q)}, and their catastrophe map, %(M), 
into the control planes, C — {(m 2 ,Q)}, for V3 Model. The dash-dotted lines in M denote 
stability-change points E, and the dash-dotted lines in C denote their catastrophe maps x(M). 
Solid lines in M (on the light-cyan colored surface) and dashed lines (on the light-magenta 
colored surface) represent stable and unstable solutions, respectively. The arrows indicated by 
"thin" and "thick" show the thin-wall limit, uj 2 — > Cj 2 ^^ = fh 2 — 1/2, and the thick- wall limit, 
ui 2 — > rh 2 , respectively. In the regions denoted by S, SU and N on C, there are one stable 
solution, one stable and one unstable solutions, and no equilibrium solution, respectively, for 
fixed (m 2 , Q). 

As m 2 becomes smaller, Qmax and Q m \ n come close to each other, and 
finally merge at m 2 ~ 0.26, below which there is no stable solution. 

The above results for the two models are consistent with the previous results for 
some special cases such as the thin- wall limit, the thick- wall limit and bubbles with 
Q = 0. 
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Fig. 2. The same as Fig. [T] but for T4 Model. Because the structure of M is complicated in this 
case, we show two pictures of M: The left one shows the upper (front) sheet of the equilibrium 
space, while the right one the lower (back) sheet. In the regions denoted by N, U, SU and 
SUU on C, there are no equilibrium solution, one unstable solution, one stable and one unstable 
solutions, and one stable and two unstable solutions, respectively, for fixed (m 2 ,Q). 

Although we have investigated only two concrete models, taking account of the 
fact that the stability structure falls into two classes in the thick- wall limit, that 
is, the fact that Q-balls are stable if n < 10/3 for the potential ([8]) and unstable 
otherwise, one expects that there are essentially two distinct stability structures in 
the general case. Then the two types of models investigated here, V3 and V4, may 
be regarded as the representatives of these two distinct stability structures. 

For example, in the gravity-mediated supersymmetry breaking model, n ) the 
lowest-order negative term of the potential is ~ — 4> 2 log 4>. Because this term cor- 
responds to n < 3 in (jHJ), the stability structure of this model falls into V3 Type. 
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Fig. 3. A schematic picture of the potential E(uj;fh 2 ,Q) of Vz Model with m 2 = 0.2 near the 
equilibrium solutions, and the locus of equilibrium solutions on (Q,E) plane. The solid and 
dashed lines represent stable and unstable solutions, respectively. 

Furthermore, because the potential is positive everywhere, which corresponds to 
rh 2 > 1/2 in Fig. 1, all equilibrium solutions are stable in this model. 

In summary, we have proposed a new method for analyzing the stability of 
Q-balls using catastrophe theory. An essential point is that, although the Q-ball 
system ([2]) includes infinite degrees of freedom, practically it can be regarded as 
a mechanical system with one variable, iv, near equilibrium solutions. Therefore, 
we have applied catastrophe theory, which was established for mechanical systems 
with finite degrees of freedom, to the Q-ball system. A similar analysis but on 
the stability of exotic black holes was done by Maeda et al. 12 ^ some time ago, and 
catastrophe theory was found to be very useful. Thus it seems worthwhile to consider 
the application of catastrophe theory to other cosmological (gravitating) solitons such 
as gravitating Q-balls, 13 ^ topological defects, and branes. It may be also interesting 
to apply the catastrophe-theoretic approach to non-relativistic atomic Bose-Einstein 
condensates, 14 ^ where Q-ball-like solitons appear. 

We thank H. Kodama, K. Maeda, K. Nakao, V. Rubakov, H. Shinkai, T. Tanaka 
and S. Yoshida for useful discussions. A part of this work was done while NS was vis- 
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Fig. 4. The same as Fig. El but for V 4 Model with m 2 = 0.3. 
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Appendix A 

Basic Idea of Catastrophe Theory 



To illustrate the basic idea of catastrophe theory, 15 ) we consider a system with 
one behavior variable x, two control parameters p,q and a potential F(x;p,q). An 
equilibrium point of x is determined by dF/dx = for each pair of (p, q). The set of 
the control parameters, C = (p,q), spans a plane called the control plane, and the 
set of equilibrium points, 



M 



(x,P,q)\f(x,p,q) = 



dF 

dx 



(A-l) 
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is called the equilibrium space. Because equilibrium points are stable if df /dx > 0, 
the boundary of stable and unstable equilibrium points are given by the curve, 

S=Ux,p,q)\^=0, f = o\. (A-2) 

The catastrophe map is defined as 

X- M^C, (x,p,q) -> (p,q). (A-3) 

According to Thorn's theorem, depending on the configurations of the image x(^)> 
all mechanical systems with stability-change are classified into several catastrophe 
types. If the number of control parameters is two, as is this example, possible 
catastrophe types are fold catastrophe and cusp catastrophe. As we show in the text, 
Q-ball models are also classified into these two types. 

If the potential F(x;p, q) is known, it is easy to find equilibrium points and their 
stability. However, even if we do not know the explicit form of F(x;p,q), we can still 
find E by analyzing equilibrium points as follows. The Taylor expansion of f(x,p, q) 
in the vicinity of a certain point P(xo,po, go) in M, where / = 0, up to the first order 
yields 

^.^-(l)-'^-^!^, if |^o. (A.4) 

Because df /dx = in U, it follows from (|A-4|) that dq/dx = in S. Similarly, 
unless df /dx = 0, dp/dx = in U. Therefore, surveying the points with dp/dx = 
or dq/dx = in the equilibrium space M, we can obtain the set of stability-change 
points E. 
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